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Abstract 

For the low-energy anti-neutrino reaction, 9 e +p —> e + + n, which is of great 
current interest in connection with on-going high-precision neutrino-oscillation 
experiments, we calculate the differential cross section in a model-independent 
effective field theory (EFT), taking into account radiative corrections of order a. 
In EFT, the short-distance radiative corrections are subsumed into well-defined 
low-energy constants the values of which can in principle be determined from 
the available neutron beta-decay data. In our low-energy EFT, the order-a 
radiative corrections are considered to be of the same order as the nucleon 
recoil corrections, which include the "weak magnetism" contribution. These 
recoil corrections have been evaluated as well. We emphasize that EFT allows 
for a systematic evaluation of higher order corrections, providing estimates of 
theoretical uncertainties in our results. 



1 Introduction 



Low-energy anti-neutrinos from nuclear reactors are well suited to determine the 
neutrino mixing angle #13, which is important for the search of CP violation in the 
leptonic sector; see, e.g., Refs. [HE]- The Double-Chooz [3], Daya Bay jl], and RENO 
Collaborations [5] are aiming to measure #13 with very high precision with the use of 
z/ e 's produced in nuclear reactors. The present upper bound to this quantity reported 
by the Chooz [3] and MINOS [6] Collaborations is: 13 < 11.4°. 

The Double-Chooz, Daya Bay and RENO experiments monitor the inverse beta- 
decay reaction on a hydrogen target 

^ e + V — > e + + n (1) 

for a known anti-neutrino energy flux. The positron yield is measured as a function 
of the positron energy. An accurate extraction of the mixing angle 613 from an 
analysis of the measured positron yield requires a precise knowledge of the radiative 
corrections (RCs). In earlier papers [3 El [9], the relevant RCs were evaluated in the 
theoretical framework developed by Sirlin and Marciano [TOl ITT] . In this framework, 
to be referred to as the S-M approach, the RCs of order a are decomposed into 
so-called outer and inner corrections. The outer correction is a universal function 
of the lepton energy and is independent of the details of hadron physics, whereas 
the inner correction is influenced by short-distance physics and the hadron structure. 
The inner corrections coming from 7 and weak-boson loop diagrams are divided into 
high-momentum and low-momentum parts. The former is evaluated in the current- 
quark picture, while the latter is computed with the use of the phenomenological 
electroweak-interaction form factors of the nucleon. Although the estimates of inner 
corrections in the S-M formalism are considered to be reliable to the level of accuracy 
quoted in the literature, the possibility that these estimates may involve some degree 
of model dependence is not totally excluded. 

We present here a calculation of the RCs to order a based on effective field the- 
ory (EFT). We use heavy-baryon chiral perturbation theory (HBxPT), which is an 
effective low-energy theory of QCD, see e.g. Ref [12] . In HB;\PT the short distance 
hadronic and electroweak processes are subsumed into a well-defined set of low-energy 
constants (LECs). In other words, these LECs systematically parameterize the inner 
corrections of the S-M approach. Therefore, insofar as there are enough sources of in- 
formation to determine the values of these LECs, HBxPT leads to model-independent 
results with systematic estimates of higher-order corrections. The use of HB%PT to 
calculate electroweak transition amplitudes for the nucleon and few-nucleon systems 
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were pioneered in Refs. [121 [EH [32], an d subsequently there have been many impor- 
tant developments. In Ref. [TB], we presented the first ever EFT-based calculation 
of RCs for the neutron /3-decay process, n — > p + e~ +v e . Because in HB^PT the 
nucleons are treated as point-like, it is expected on general grounds that the order-a 
RCs are common between neutron /3-decay and inverse /3-decay. Meanwhile, it should 
be mentioned that, in the counting scheme adopted here and in Ref. [16]. the order-a 
RCs are of the same order as the m" 1 nucleon- recoil corrections including the "weak 
magnetism" contributions^] and hence a consistent EFT calculation should include 
these recoil corrections simultaneously. We present here such an EFT calculation, 
taking advantage of the fact that the m~ 1 -expansion is a natural part of our counting 
scheme and thus dictates how to incorporate recoil corrections order by order (see 
later in the text). 

Since exactly the same LECs are involved in the EFT calculations of inverse /3- 
decay and neutron /3-decay, we can in principle use the existing neutron /3-decay 
data to determine those LECs and make a model-independent estimation of RCs for 
the inverse /3-decay, provided that the m" 1 recoil corrections are properly taken into 
account. In this connection, we note that an attempt has been made in the litera- 
ture [19] to directly relate the neutron decay rate with the inverse /3-decay cross 
section, assuming that the order-a corrections (RCs and recoil corrections combined) 
are common between these processes. As mentioned, this assumption is justified as 
far as the genuine RCs of order a is concerned. However, as described later in the 
text, our calculation shows differences between the m" 1 corrections for inverse 0- 
decay, Eq. ([T]), and those for neutron /3-decay. We therefore caution against writing 
the cross section for the reaction in Eq. (TjQ) in terms of the neutron mean life, r n , as 
advocated in Refs. f7\ [T9]. 

This paper is organized as follows. In section 2 we explain a theoretical framework 
to be used and present the results for the order-a RCs. In section 3 we consider the 
recoil corrections and compare our results with an earlier work[19j. Section 4 gives a 
summary of our calculations and conclusions. The appendix describes some technical 
details concerning the HB;\RT treatment of the infrared singularity. 

2 The QED corrections 

We use here essentially the same theoretical framework as in Ref. [16], in which we 
calculated RCs for neutron /3-decay. We therefore give only a brief recapitulation of 

1 The importance of the nucleon-recoil corrections was emphasized by, e.g., Vogel and Bea- 
com [19]. 
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our formalism, relegating details to Ref. 

Our calculation is based on the (Q/A x -expansion scheme, where Q 



E u —A N —m e 



(An = m n —m p ) represents a typical four-momentum transfer for incident low-energy 
reactor anti-neutrinos (E v < 10 MeV), and A x ~ Anf^ « 1 GeV (/ w = 92.4 MeV 
is the pion-decay constant) is the chiral scale. It is to be noted that the expansion 
parameter in our scheme is very small and that, as explained in more detail below, 
the lowest order recoil corrections ~ Qj m N are of the same order as the lowest order 
radiative corrections; viz., Q/m N ~ a/(27r) ~ <3/A x ~ 1CT 3 , where m N = (m p +m n )/2. 

The leading order (LO) transition matrix element for the inverse /3-decay, Eq.fll]), 
is evaluated ignoring nucleon recoil and radiative corrections. The next-to-leading 
order (NLO) corrections in our counting scheme are the recoil corrections (~ Q/m N ) 
and the radiative corrections (~ a/(27r)). The recoil corrections, which include the 
"weak magnetism" term, will be specified in Eq.(j5]) below. For the sake of the trans- 
parency of presentation, we shall in this paper separate these corrections from the 
m^ 1 (kinematic) corrections to the phase-space. 

The effective lagrangian relevant to our calculation includes the relativistic lep- 
tonic weak interaction current and the LO and NLO heavy-baryon lagrangian 
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Cqed m Eq.Q is the usual QED lagrangian, where F^ v 
d^ + ieA^ is the covariant derivative; for the gauge parameter £ A , we use here £4 = 
1 (Feynman gauge). Cnn is the heavy- nucleon lagrangian including the photon- 
nucleon interaction, and Cnn4>4> is the low-energy LO and NLO current-current weak 
interaction. We give in Eq.(jSJ) the explicit forms of NLO nucleon-recoil terms dictated 
by HBxPT. In the above, g A = 1.267 is the axial coupling constant, while is the 
nucleon velocity vector, and is the nucleon spin; they satisfy v ■ S — 0. We choose 
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here = (1,0) and = (0, a/2). In the NLO part of the lagrangian the nucleon 
isovector magnetic moment /i v = fi p —fi n = 4.706. The low-energy constants (LECs), 
ei, e 2 , e v and counter-terms which regulate the ultraviolet (UV) divergences of 

the virtual photon-loop diagrams. These LECs incorporate the short-range radiative 
physics that is not probed in a low-energy process. The LECs, e\ and e2, are related 
to the wave-function renormalization factors of the positron and proton, respectively. 
The LECs, e v and e A , are related to the Fermi and Gamow- Teller amplitudes. The 
Fermi coupling constant, Gf = 1.166 x 10~ 5 GeV -2 , is determined from muon decay, 
and the CKM matrix element, \V ud \ = 0.97418 ± 0.00027, is given by the PDG [20]. 

For later convenience, we rearrange the LECs in Eq. (IH]) by rewriting the hadronic 
part in the first line in Eq.([5]) as 

NT + [v»-2g A S> x ]N + (^je Y NT + [v»-2~g A S»}N + 0(a 2 ) , 

where we have introduced the redefined axial coupling constant, g A = g A [l+J^(e A — e v )]. 
As in the neutron /3-decay case [IE] , to the order of our concern, g A can always be 
replaced by g A . This also applies to the NLO recoil contributions since the m" 1 recoil 
corrections are of the same order as order-a corrections in the adopted counting 
scheme. The order-a radiative corrections to the nucleon magnetic moments are for 
the same reason of higher order in our scheme and hence neglected in this work. 

In this paper we derive a model-independent expression for the lowest order ra- 
diative and recoil corrections to the reaction, u e (Pu) + p(Pp) ~~ e+ (Pe) + n (p n ), where 
the four-momentum of each particle is indicated in the parentheses. We shall concen- 
trate on an experimental setup in which none of the particle spins are monitored by 
the detector. There is one subtle aspect of the above reaction which deserves some 
discussion. In experiments, the final state positron will always be accompanied by 
(often undetected) soft bremsstrahlung photons. If the bremsstrahlung photon en- 
ergy, E y is less than the detector resolution, A, the energy recorded by the detector is 
the sum of the actual outgoing positron energy, E e , and the bremsstrahlung photon 
energy, £7 7 ; i.e., E = E e + J5 7 is what is measured as the "detected positron" energy, 
with the corresponding "detected positron" momentum being \p\ = \J E 2 — m*. The 
two processes we evaluate are u e + p — > e + + n and v e + p — > e + + 7 + n. Due to 
the finite detector resolution the second bremsstrahlung process is not observed; it 
only contributes to the RCs of the first process, i.e., the soft bremsstrahlung photons 
are an integral part of the "detected positron". Thus, the first process has become, 
Ve(Pv) + p{Pp) e + (p) + n(p n ), where the positron momentum p e has been replaced 
by p in order to indicate that the soft bremsstrahlung process has been incorporated 
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into this "effective" reaction. The cross section for this "effective" reaction is given 
in terms of the effective invariant amplitude M.: 

f(E)\(l + 3g 2 A )g 1 (f3) + (l-g 2 A )g 2 (f3)f3cos9 e \d(cos9 e ), (6) 



V V2 

where /3 = \p\/E = J ' E 2 — m\jE is the velocity of the outgoing "detected positron" 
for a given incident (anti-)neutrino beam energy, E v and detector reading, E; cos(8 e ) = 
p v - p, and f(E) is the phase-space factor to be discussed later in the text (see, 
Eq. (TT7|) ). The two velocity-dependent functions, Gi(f3) {i = 1, 2), are written up to 
NLO as 

gtf) = 1 + sr d (P) + — Q? cm \[i) , (7) 

Here G™ d (f3) (see, Eqs. ( 1T0|) and f|TT]) ) represent the lowest-order radiative corrections, 
and Q^ ecod ((3) (see, Eq. f|T8|) ). which will be evaluated in the next section, represent 
the recoil corrections arising from the lagrangian in Eq.()5]). The calculation of the 
function §l ad (f3) is described next. 

For the analysis of the radiative corrections, we explicitly distinguish between the 
outgoing positron and the bremsstrahlung photon. There are two distinct categories 
of radiative corrections, the bremsstrahlung and the virtual photon loop corrections. 
The corresponding Feynman diagrams are shown in Fig. [TJ Since O(a) and 0(m~^) 
are of the same order in our counting scheme, the invariant matrix element, Aibr for 
bremsstrahlung is evaluated assuming E n = m n . The differential cross section for the 
radiative process, i'e(Pu)+p(Pp)^e + (Pe)+n( y p n )+'j(kj^, is given by 

da br (u e p ->■ e + nj) 



d 3 p e d 3 k 
8m p m n E u J (2vr) 3 2£ e (27r) 3 2£, 



7 

x (2ir)S(E v -A n -Et- E e ) 1 - £ \M br \ 2 , (8) 

spins 

where E = E e + E 1 = E v — A^v is the maximum energy of the "detected positron" 
in the static nucleon approximation, i.e., m e < E < E. The bremsstrahlung matrix 



The four- moment a, p e and k denote momenta in the static nucleon approximation, i.e., p e = 
0{m~^ 1 ), k = k — ©(m^ 1 ), and correspondingly, E = E — ©(m^ 1 ). 
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element squared with the static neutron is 
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E e E u (p e -k), 
m\E^ \ 
E e E u {p e -~k)) 



(9) 



We note that the above expression for \Mbr\ 2 is identical to that for neutron {3- 
decay derived in Ref. [33]. We also remark that Eq.([9]) was derived earlier by Fukugita 
and Kubota [9] , who used the S-M approach [TQl E] and a finite photon mass in order 
to regulate the infrared (IR) singularity. In the integration over the bremsstrahlung 
photon energy in Eq.® the maximum photon energy occurs when E e = m e , i.e., 
E™ ax = E — m e Jf| In this context, the same question again arises as to whether 
or not the experiment can distinguish between the two final states, n + e + + 7 and 
n + e + . If the detector resolution in the experimental setup is such that one can detect 
photons with an energy in the interval A < < E™ ax , we should integrate the 
bremsstrahlung photon energy E y = \k\ from to A in Eq. ([8]). However, if the 
experiment is unable to distinguish these two final states, we should integrate from 
to E™ ax . In order to compare with earlier works, we concentrate here on the latter 
case. The integral over the radiative photon spectrum invariably gives rise to an IR 
singularity. We use dimensional regularization to deal with the IR singularity; some 
details regarding the bremsstrahlung integral are presented in the appendix. As is 
well known, the IR singularity appearing in Eq.® is cancelled by the contributions 
from virtual photon-loop diagrams in accordance with Bloch and Nordsieck [17], see 
also [18J. The evaluation of the loop diagrams in dimensional regularization can be 
found in the literature, see, e.g., Ref. [IE]. It is notable that, apart from the so-called 
"Coulomb factor" it 2 / (3, which arises in, e.g. neutron /3-decay from a photon-loop 



3 The approximate integrals considered in Ref. [7] give the analytic expression in Ref. [5] provided 
the lower limits of the integrals are changed from 1 MeV to m e . 
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diagram, the matrix element given by these virtual photon loops is identical to the 
one in neutron /3-decay. 

The UV-divergencies originating from the photon loop diagrams are regulated by 
the LECs in the lagrangian. These LECs are renormalized by the usual effective field 
theoretical method based on dimensional regularization of the loop integrals, see e.g., 
Refs. [T2"} [22] • The finite LECs renormalized at the scale \x are 

,2 \ 



) = ey,A 



1, .3 
2 (6l + 62) + 2 



d-4 



7s + ln(47r) + 1 




The LEC, (/i 2 ), which was introduced by Ando et al. [15] in the evaluation of the 
RC for neutron /3-decay, subsumes short distance physics not probed at low energies 
and depends on the regularization scale \x. 

Combining the bremsstrahlung and virtual photon-loops contributions calculated 
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to order a, and noting that j3 - 
0(m- N l ) contributions, Q\ ad {fi) 
terms of 0[a.m~^\ we choose to write the results in the following form: 
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where the "inner" corrections, which are independent of /3, are encoded in the LEC 
ey(/x 2 ) and defined as ey(/z 2 ) = (/z 2 ) + |. The "outer" radiative corrections consti- 
tute the well-known, model- independent, long-distance QED corrections that do not 
contain any hadronic effects, and are given by 



<W/3) 
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The above expressions for 5 out and S out reproduce the results obtained by Fukugita 
and Kubota [9J. We also note that 5 out = h(E,E ), where h(E,E ) is the function 
introduced by Sirlin [23] . 

As mentioned, e^(/i 2 ) also appears in the expression for the RC for neutron /3- 
decay. Therefore, it is in principle possible to determine e^(fi 2 ) using relevant high- 
precision low-energy data involving baryons. Due to lack of useful experimental data, 
Ando et al. [16] determined (/i 2 ) at n — m N by comparing their results for neutron 
/3-decay with those obtained in the S-M approach [TUJ [TT] . 



3 The m recoil corrections 

N 

As mentioned, these corrections have two different origins. One comes from the 
lagrangian itself, and the other arises from the expansion of the kinematic factors 
in the phase-space integral. Below we treat these two types of recoil corrections 
separately and compare our results with those in Ref. [19]. It is to be noted again 
that, in evaluating the 0{mjl) corrections, we can neglect O(a) radiative effects, 
since 0(a m^ 1 ) terms are of higher order in our counting scheme. One can, therefore, 
assume that the outgoing positron energy, E e « E, and correspondingly, the positron 
velocity, /3 e ~ /3. 



Kinematic (phase space) corrections 



The phase space factor, f(E), appearing in Eq.(JBJ to the lowest order (LO) in the 
m^ 1 expansion is given by f(E) = E 2 (3/tt with the neutron regarded as being static, 
i.e., E n = m n + (9(m~ 1 ). To NLO, the above expression for f(E) needs to be corrected 
to incorporate the kinetic energy of the recoil neutron from the relation E n = m n + 
(p u — p e ) 2 / (2m N ) + 0(m~ 2 ). Corresponding to this change in E n , we have 



E 



E u -A N -(p u -p e ) 2 /(2m n ) + 
E 



1 / ~ A 2 — 

1 EJl-(3cos9 e ) + 

m N \ 2E 
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where, as earlier, E — E v — Ajv, and the positron velocity becomes 
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where (3 = y E 2 — ml/ E. Note that the positron energy, E and the velocity, j3, are 
equal to the recoil-corrected E^ and in Ref. [19], respectively. Reflecting these 
changes, the phase space integral in Eq.flSD needs to be corrected as follows: 



f(E) fm n + E n 



An 



E n 



6{7) 



dE 



;i + 3^) + (l-^)/3cos^ e , (16) 



where F = E u — A N — E — (p u — p e ) 2 / (2m N ) H . The factor (m n + E n )/E n ~ 2 in 

Eq. (IT6"|) has corrections of order m~ 2 , and the Jacobian factor produces the following 
NLO phase space factor in Eq.Q: 



f(E) 



E 2 (3 



71 
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(17) 



where the expressions for E and (3 are given in Eqs. (IT3|) and ( TT5]) . 



Corrections from the next-to-leading-order lagrangian 

The m" 1 corrections to the Lagrangian are explicitly written in Eq.(JS]). As noted 
before, the radiative corrections to these additional terms in the Lagrangian are of 
higher order than NLO in our counting and hence need not be considered here. Eval- 
uating the NLO lagrangian recoil correction contributions, illustrated in diagram (a) 
in Fig. [U we obtain the recoil terms in Eq. (JTj) 

grecoii^ = ^2 e ( 1 - 2 9a^ + gl \ £ ( 1 + 2g A fi v + g\ \ 



g™*{p) = e 1 " A ^ v 9 1 - e, 



■2\ /i ■■),,,, ,:l 

2 I " v \ 1 



,recoil fQ \ _ tpI 1 + 2 9^v + 9a\ p ( 1 ~ 2 9a^v + 91 



i-gi J V 1-5 



A 



Comparing these results with those obtained for neutron /3-decay [IB], we note that 
there are several relative sign differences^ Apart from the m^ 1 phase-space cor- 
rections in neutron /3-decay, the mj^ 1 corrections (arising from the Lagrangian m^ 1 
interaction terms) relevant to the neutron life-time are contained in the Cq factor in 
Eq.(14) of Ref. [16]. Noting that in neutron beta-decay E™ ax = E u + E e + O^ 1 ), 
we may rewrite the C factor as 

'l + 2fi v g A + g 2 A \ , „ fl-2fi v g A 



C (E e ) = 1 + IpEA 7^ 2 ' 1 4- E v i 

m N { \ l + 3g 2 A J \ l + 3g 2 A 




4 This is in contrast to the order-a RCs which are universal at NLO in effective field theory. 
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Comparison of Cq with Q[ ecml ((3) in Eq.f llSp clearly indicates that the m" 1 recoil cor- 
rections are not identical for the neutron /3-decay and the inverse /3-decay. Moreover, 
since the weak- magnet ism term is dominant, the difference between Cq and G[ ecml (f3) 
are of the same magnitude as the corrections themselves. 

Combining the m" 1 Jacobian factor in the square brackets in Eq. f lTT]) . and the re- 
coil correction arising from the lagrangian, Eq. ( TTgj) . we confirm the recoil corrections 
given in Eqs. (12) and (13) in Ref. |19] . We prefer to keep these two m' 1 corrections 
separate since one is of a kinematical origin (phase space correction), whereas the 
other is of a dynamical origin arising from the transition matrix element. 

4 Discussion 

In this paper we have derived a model-independent expression for the radiative and 
m^ 1 corrections for the low-energy anti-neutrino proton reaction to next-to-leading- 
order in an effective field theory approach. We have shown that short-distance physics 
not probed in this low-energy reaction can be subsumed into a single low-energy con- 
stant e^(/i 2 ). In the QM x -expansion scheme adopted here, the 0(a) and C^Q/m^) 
corrections are considered to be of the same order for the reactor anti-neutrino en- 
ergy range. We have found that the m~ l corrections appearing in Eq. ( TT8|) . which 
originate from the lagrangian Eq. fl5]), are different from the m~ l corrections found 
in neutron /3-decay, see e.g., Ref. [IB]. Therefore, to the order under consideration, it 
is not advisable to write the inverse beta-decay cross section (or the positron yield) 
in terms of the neutron mean life r n , as advocated in Ref. [TJ. 

The short-distance hadronic physics associated with the LEC, (/i 2 ), was exten- 
sively discussed in Refs. [101 HI] - The processes involved in ey(/i 2 ) were studied from 
an effective field theoretical perspective in Ref. [32]. In principle, we should be able 
to determine the LEC, e^lfi 2 ), from high-precision experimental data. Relegating 
this determination to future study, we choose here to estimate (/x 2 ) at the scale 
\x = m N by comparing the short-distance radiative corrections calculated in the S-M 
approach [10J [TTJ and the expressions obtained in EFT in Ref. [16]. The result is 

e£(m 2 ) = A\n(^j+\n(^j+2C + A g = 18.31 - .25 + 1.78 - 0.34 = 19.5, (19) 

where, for the sake of definiteness, the value of the axial matching mass 771^=1.2 GeV 
has been used although its value involves uncertainty |T0l [11]. With this value of 
ey(m 2 ), the correction term involving LEC in Eqs. f lTOj) and ffTTl) is estimated to be 
(a/27r) ey(m 2 ) ~ 0.023. The dominant first term in Eq. ( 1191) arises from well-known 
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additional box diagrams with Z-exchange, replacing the photon-exchange, in electro- 
weak theory (TUJ E]. This electro- weak physics can be naturally included in our 
approach. However, for an easy comparison with the neutron beta-decay radiative 
corrections evaluated in Ref. [16], we prefer to keep this contribution in the above 
LEC. As for the last two terms in Eq. f fT9|) . we remark that A g involves genuine 
short-distance hadron-structure physics, whereas the constant C arises from photon- 
loop diagrams in which the photon couples to the nucleon magnetic moments and 
also from the hadronic form factors. The long-range parts of these corrections are 
naturally included in EFT at higher orders than considered in this paper. 

As a final comment we note that in our work we have used the value of the 
Fermi constant Gp determined from the muon lifetime measurement. The theoretical 
expression for Gp is evaluated in standard electroweak theory, and it naturally in- 
cludes log-terms involving mz- These log-terms appear in our expression for ey(/i 2 ), 
Eq. ( 119]) . and was also considered in Ref. [16], see e.g., Refs. [TTJl [TTJ for details. 

In summary the integrated cross section for reaction ([1]) is 

o = (G F V ud ) 2 ^(l + 3g 2 A )(l + ^-gi ad 0^ 



where as before, E = E v — (m n — m p ) and ft = yE — m^fE, and where all m^ 1 
corrections in Eq. (l20|) except Q[ ecml (/3) of Eq. (fT8|) originate from the phase-space 
factor f(E), Eq.flTT]). 
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Appendix 

We use dimensional regularization to isolate the IR singularity. The 3-dimensional 
integral over k in Eq.® is replaced with a d = 4 — 2e dimensional integral where 
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e < for the purpose of handling the IR singularity, i.e. Eq.flS]) is rewritten as 

- 1 J^—L r m d\~k\ \~k\ d - 2 d d - 2 n~ AM - V \M h I 2 (2D 



where p e = y (E— \k\) 2 — m 2 . We note that in dimensional regularization, the angular 
integration / d d ~ 2 f2^ yields (p e -fc = cos 9k) 



^A-d ,d d ~ 2 fi^[l-cos 



(2tt)^ J [1 - /3cos(^,)] 2 

,4-d ( 

\ [l + |e| (7 B -ln(47r)) 



. £ 2 /? 3 U-/3, 



+ CC9) + 0(e 2 ) 
P 2 



where (3 = \J E 2 — m 2 /E and the function C(/3) is given by (see e.g. Refs. [2T] ) 

+ 21n2 



i + im(I±£ 

2/3 V1-/3, 
1 



1-Iln(i±£ 



2(3 \1-P, 



2/3 



in 



(22) 



and L(x) is the Spence function 

L(x) = —Lii(x) = / 
Jo 



* , ln(l — t) 



/o t 

The integral over the photon momentum / dk k d ~ 5 oc l/|e| exhibit the IR-singularity. 
When we combine our expression for the integrated bremsstrahlung cross section with 
the contributions from the virtual photon loops, we find that the IR singularity is 
removed as it should. 
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Figure 1: Feynman diagrams contributing to the cross section which includes 0(a) 
QED corrections: (a) - leading order (LO) Born amplitude; (b) and (c) 
bremsstrahlung amplitudes; (d), (e) and (f) - - virtual photon-loop diagrams; (g) 
and (h) — counter-term amplitudes involving the LECs e\ and e2, respectively; (i) - 
counter-term amplitudes involving both of the LECs, e v and e A . The m" 1 correction 
is represented by diagram (a) but with the use of the vertex that arises from the 
m^ 1 (NLO) part of the lagrangian in Eq.(jSJ). To the order of our concern we do not 
consider diagrams in which a photon couples to the nucleon magnetic moments since 
this is 0{a/m N ) or NNLO. 
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